Nodes in the Gap Function of LaFePO, the Gap Function of the Fe(Se,Te) Systems, 

and the STM Signature of the s± Pairing 
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We reiterate, in more details, our previous proposal of using quasi-particle interference to de- 
termine the pairing form factor in iron-based superconductors. We also present our functional 
renormalization group(FRG) results on LaFePO and Fe(Se,Te) superconductors. In particular we 
found that the leading pairing channel in LaFePO is nodal s±, with nodes on electron Fermi surfaces. 
For Fe(Se,Te) system we found fully gapped s± pairing, with substantial gap anisotropy on electron 
Fermi surfaces, and large gap is concentrated in regions with dominant xy orbital character. We 
further fit the form factor obtained by FRG to real space orbital basis pairing picture, which shows 
more clearly the differences between different iron-based superconductors. 

PACS numbers: 74.20. Rp,74.20.Mn 



In the midst of many remaining issues of the iron-based 
superconductors, the "pairing symmetry" has attracted 
a lot of attention. From the theoretical side, the answer 
converged rather quickly to the s± form factor 1-7 . How- 
ever, on the experimental side conflicting evidences exist 
for the presence of nodes 8 and full gap 9 . Note that we 
use "form factor" rather than "pairing-symmetry" to de- 
scribe s±. This is because from the symmetry point of 
view there is no difference between the s± and the usual 
s-wave pairing. Indeed, under the action of the crystal 
point group, both transform as the identity (trivial) rep- 
resentation. This is in marked contrast with the d x 2_ y 2 
pairing symmetry of the cuprates, which transforms as a 
distinct irreducible representation upon the point group 
operations. Thus while there is sharp (symmetry) dis- 
tinction between the d x 2_ y 2 and the usual s-wave pairing, 
such distinction does not exist between s± and s. This 
is why phase sensitive measurements probing the relative 
phase of the superconducting order parameters residing 
in regions with different spatial orientation are ideal to 
rule in or rule out the d x 2_ y 2 pairing 10 ' 11 , while they can 
not definitively prove or disprove the s± form factor 12 . 

Nonetheless at this point there exist two types of pro- 
posed experiment which can address whether the gap 
function on the electron and hole pockets are indeed of 
opposite sign. Both of them provide information about 
the relative (spinor) phase of the quasiparticle wave func- 
tions on the electron and hole pockets. One of them is 
the detection of neutron resonance mode in the supercon- 
ducting state. As suggested in Ref. 13,14 , if the pairing 
form factor is s± one expects neutron resonance peaks 
at momenta (±7r, 0), (0,±7r) to occur in the supercon- 
ducting state. The energy location of the peaks is less 
than the sum of the minimum gaps on the electron and 
hole pockets. Interestingly the resonance peak has been 
observed 15-19 . 

In a recent paper 20 the present authors proposed 
the second type of experiment - STM "quasiparticle- 



intcfercncc" experiment. This experiment can also probe 
the relative (spinor) phase between the quasiparticles on 
the electron and hole pockets. Here we repeat the argu- 
ment provided in Ref. 20 . If the electron and hole Fermi 
surfaces have out-of-phase order parameters, the associ- 
ated quasiparticle Nambu spinor will be orthogonal (left 
panel, Fig. 1). 

For example, under the gauge where the order param- 
eter is real, one of them will be ~ (1,1), and the other 
~ (1,-1). The matrix elements of impurity- induced 
quasiparticle scattering from the hole to the electron 
Fermi surfaces (or vice versa) is non-zero/zero due to 
a scalar/magnetic impurity (which acts as the Pauli ma- 
trix u z /l2x2 in the Nambu space). The reverse is true 
for the scattering between two hole or two electron Fermi 
surfaces. As a result, if the main scattering source are 
scalar impurities, the quasiparticle interference peaks as- 
sociated with (±7r, ±7r) (which arise from the scatter- 
ing between two electron or two hole pockets, note that 
throughout this paper the "unfolded" zone notation of 
momentum space will be used) will be suppressed for bias 
close to the superconducting gap. In contrast, for mag- 
netic impurity scattering, the quasiparticle interference 
peak around (±7r, 0) and (0, ±7r) (which correspond to 
the scattering between an electron and a hole pocket) 
will be suppressed. The contrast between these two 
sets of quasiparticle interference peaks should diminish 
as the bias is increased/decreased from the gap because 
the spinor phase of the quasiparticle wave functions are 
no longer solely determined by the gap parameter. The 
quasiparticle interference idea 21 has been developed into 
a fruitful spectroscopy of studying the cuprates 22 . It has 
also been used method to infer the pairing sym- 

metry by Hanaguri et at 23 . In a recent beautiful STM 
experiment 24 done on the Fei+z(Se,Te) compounds, the 
authors control the relative degree between scalar and 
magnetic impurity scattering by changing the density of 
superconducting vortices (which act as magnetic scat- 
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FIG. 1: (Color online) (a) Schematic representation of the s± pairing form factor, and the associated quasiparticle Nambu 
spinors. Solid black circles represent the normal state Fermi surfaces. Thickness of colored region around Fermi surfaces 
indicate magnitude of the gap |A(k)|. Blue means positive and red means negative. The two component column vectors are 
the quasiparticle spinor wave functions associated with different Fermi surfaces, (b) A schematic representation of the pairing 
form factor for LaFePO obtained by FRG. There are four nodes of gap around each electron Fermi surface. 



tcring centers). Interestingly the above behavior is ob- 
served. 

Combining the neutron and the STM experiment it 
is fair to say that the evidences for the s± pairing is 
strong. The purpose of the present paper is to tie up 
the loose ends on the theoretical side. We study the 
pairing form factors for LaFePO and Fe(Sc,Te) systems. 
These systems exhibit superconductivity without charge 
doping, and are presumably less prone to disorder ef- 
fects. LaFePO, among all iron-based superconductors, 
shows the strongest evidence for the existence of gap 
nodes 25,26 . For FeSe, there have been contradicting re- 
ports of nodeless 27 and nodal 28 superconductivities, as 
in many of its pnictides relatives. Bulk FcTc is not su- 
perconducting but becomes superconductor under tensile 
stress as thin films 29 . 

We apply the functional renormalization group (FRG) 
method to study the pairing form factor. In particular we 
incorporate the realistic band structures of these materi- 
als in the form of tight binding models. We then trans- 
form the real space on-site interaction into two-particle 
scattering vertex functions defined in the band cigenfunc- 
tion basis. It turns out that much of the differences be- 
tween the pairing form factors of LaFePO, Fe(Se,Te) and 
the LaFeAsO systems originate from the differences in 
the Fermi surface topology as well as the orbital contents 
of the band wave functions. Thus the realistic band struc- 
ture is indispensable for our purposes. 

As pointed out in Ref. 30 ' 31 , the scattering vertex func- 
tion acquires important angle dependence around the 
Fermi pockets due to the change in the orbital content 
of the band wave functions around each Fermi surfaces, 
and this can induce strong variation in the s± gap func- 
tion around the electron pockets. Recently it has also 
been shown that a particular type of angular dependent 
scattering amplitude can even lead to nodal s± form 
factor 32-34 . 



The Hamiltonian we start with is given by the sum of 
the band structure part and the interaction part H = 
Hq + Hi n t ■ Because the relevant bands are mostly iron in 
character, we follow the literature in using a tight-binding 
model consists of only the iron orbitals. Moreover we fo- 
cus on the two dimensional X-Fc-X (X=As, P, Sc) trilaycr 
in which the Fe form a square lattice. For both materials 
five-orbital tight-binding models (1) fitted to the band 
structure are used, 

r,s ij a.b 

where a, b — 1,5 label the five Wannier d-orbitals (3z 2 — 
r 2 ,xz,yz,x 2 — y 2 ,xy) of Fe. Here s =t,i labels spin, 
i,j labels the iron sites. In the above x and y refer the 
the two orthogonal Fe-Fe directions. After a staggered 
definition of the phase of the xz and yz orbitals, the tight- 
binding Hamiltonian can be made to have the translation 
symmetry consistent with one Fe per unit cell. In the rest 
of the paper we shall use this unit cell and the associated 
Brillouin zone (which is referred to as the "unfolded zone" 
in the literature). 

For iJj nt we only consider the following on-site inter- 
actions 

H\ni = o E^ E^(E] ^ ab C ias C ibs' C ibs' C ias 
i s,s' a.b 

+ ^2 Jab - (4as C ibs + h - C )(4bs> C ias> + h - C : ] ■ 
a^b 

(2) 

where :: means normal ordering, h.c. means hermitian 
conjugate of the previous term. When the parameters 
are available, we have studied the full-orbital dependent 
interaction 35 . However in the rest of the paper we shall 
focus on the simpler case where U aa = U, U a b = U' 
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for a b, and J a b = J for a ^ b. For cases where 
the orbital-dependent interaction parameters are avail- 
able we have checked that the above simplification does 
not change the qualitative nature of the results. The 
functional renormalization group procedure and result 
analysis exactly follow that used in previous papers by 
the same authors 5 ' 20,31 , except that we have achieved 
slightly better momentum resolution in the present work. 
In the following we present the results for LaFcPO, FeSe 
and FeTe separately. 

LaFePO: For LaFePO there is a discrepancy between 
the results of quantum oscillation experiment 36 and the 
band structure calculations 32 ' 37 . While the band struc- 
ture calculations predict the presence of a three dimen- 
sional 3z 2 — r 2 orbital hole-like Fermi surface centered 
around k = (tt, tt, tt), the quantum oscillation experiment 
did not observe it. In our study we follow the quantum 
oscillation experiment and adapt a band structure model 
where there is only two-dimensional-like Fermi surfaces 
associated with the two hole pockets around k = (0, 0) 
and two electron pockets around k = (tt, 0) and (0,7r), 
respectively. It is important to note that unlike the hole- 
doped 122 compounds, there is no two-dimensional xy 
orbital character hole pockets centered at k = (tt, tt) (See 
Fig. 2, left panel). We shall also comment on why we do 
not believe the existence of the (tt, tt, tt) hole pocket from 
the theoretical persepctive. 

The FRG flow of the interaction associated with pair- 
ing and the the (tt, 0) spin density wave (SDW) arc shown 
in the middle panel of Fig. 2. (Strong negative interac- 
tion implies ordering instability.) Despite the fact that 
at high energies the SDW interaction is more negative, 
the superconducting pairing overwhelms the SDW inter- 
action at low energy cutoffs. In the right panel of Fig. 2 
we present pairing form factors associated with the four 
Fermi surfaces in the left panel. We note the following 
two facts. (1) The pairing form factor on each electron 
Fermi surface has 4 nodes. This is schematically repre- 
sented in the right panel of Fig. 1(b). (2) The mean gap 
function on the electron pockets has opposite sign from 
that on the two T hole pockets. Hence it is justified to 
call them s± pairing. The presence of gap nodes on the 
electron pocket is consistent with the observation of the 
linear-T dependent penetration depth 25,26 . 

Next we address the issue of the absence/presence of 
the 3D hole pocket around (tt,tt,tt), i.e., the conflict be- 
tween the quantum oscillation experiment and the band 
structure results. A cross section of the band structure 
calculation results at k z = tt is shown in the left panel 
of Fig. 3. The Fermi surface #3 is the controversial hole 
pockets. It is primarily made up of the 3z 2 — r 2 or- 
bital. The rest of Fermi surfaces are connected to those 
at k z = without significant dispersion. Our FRG result 
shows that while the bare pair scattering between the 
electron and the 3z 2 — r 2 -like hole pocket is appreciable, 
as the energy cutoff is progressively lowered, such pair 
scattering diminishes. This is shown in the right panel 
of Fig. 3, where the black diamond symbols mark the 



form factor on the #3 Fermi surface - it is vanishingly 
small 38 . Clearly if this ungapped Fermi surface exists, it 
would change the temperature dependence of the pene- 
tration depth entirely. Thus we support the quantum os- 
cillation results and believe that the (tt, tt, tt) hole pocket 
is absent in LaFePO. This makes the band structure of 
LaFePO essentially two dimensional like the other 1111 
and 122 compounds. 

FeSe: Experimental samples of FeSe usually con- 
tain small amount of excess Fe. But they are experi- 
mentally shown as non-essential or even destructive to 
superconductivity 39 , and will be ignored in our study. 
The Fermi surface of the FeSe tight-binding model is 
shown in Fig. 4. Like the doped 122 systems there are 
two hole pockets centered around k = (0, 0), two electron 
pockets around k = (tt,0), (0,tt), and one hole pocket 
around k = (tt, tt). 

The FRG flow of the interaction associated with pair- 
ing and the (tt, 0) spin density wave (SDW) are shown in 
the middle panel of Fig. 4. Again, while the SDW inter- 
action is stronger (more negative) at high energies, it is 
surpassed by the superconducting pairing as the energy 
cutoff is lowered. In the right panel of Fig. 4 we present 
pairing form factors associated with the five Fermi sur- 
faces in Fig. 4. We note the following facts. (1) The 
Fermi surface shapes and orbital contents(not shown) are 
very similar to that of LaFeAsO. (2) The gap function on 
the two central hole pockets is quite small compared with 
that on the (tt, tt) hole pocket (mainly xy orbital) and the 
electron pockets. This suggests the main pairing source 
is the Cooper scattering between the (tt, tt) hole pocket 
and the electron pockets. (3) The gap on the electron 
pockets has substantial variation similar to our previous 
results for the 1111 systems 5 . However in this case the 
large gap is concentrated in regions with dominant xy or- 
bital content, i.e. the portion of electron Fermi surfaces 
facing the central hole pockets. This is different from our 
previous result of LaFeAsO where the xy orbital domi- 
nant part of Fermi surfaces has minimal gap. (4) The gap 
function takes on opposite sign between the electron and 
hole pockets. Hence the pairing form factor is s±. This 
is consistent with STM quasiparticle interference results 
of Hanaguri et al. 2i . 

FeTe: Bulk FeTe is not superconducting and shows a 
different antifcrromagnetic(AFM) order than the parent 
iron pnictides 40 . Its optical conductance does not show 
a clear Drude peak and does not exhibit a significant 
change across magnetic transition 41 . The magnetic mo- 
ment in AFM state is large, about two Bohr magneton 40 . 
All these facts indicate that FeTe may be very differ- 
ent from FeSe and iron pnictides, and probably more 
strongly correlated. However upon Se substitution of 
Te superconductivity appears. Interestingly, in those su- 
perconducting samples neutron scattering showns broad- 
ened (tt,0) and (0,tt) peaks similar to those of the 122 
family. This provides circumstantial evidences that the 
(tt, 0) and (0,7r) magnetic scattering is tied to the super- 
conducting pairing. Recently superconductivity in FeTe 



4 



FS 1 o FS 2 □ FS 4 a FS 5 




FIG. 2: (Color online) Left panel: the k z = Fermi surfaces of LaFePO. Middle panel: The FRG flow of the interaction 
strength associated with the (tt, 0) SDW and the superconducting pairing. A is the energy cutoff. Right panel: the pairing 
form factor. The horizontal axis 6 is the polar angle of Fermi surface points with respect to the center of each Fermi surface. 
9 = means +k x direction from the center of the Fermi surface. This notation will be used throughout this paper. 
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FIG. 3: (Color online) Left panel: the Fermi surfaces of LaFePO at k z = tt. Middle panel: The FRG flow of the interaction 
strength associated with the (tt, 0) SDW and the superconducting pairing. Right panel: the pairing form factor determined by 
FRG. 



thin films under tensile stress was also reported 29 with- 
out Se substitution. Knowing the magnetic properties of 
the superconducting Fc(Sc,Tc) systems it is reasonable to 
expect that these film to exhibit (71", 0), (0, tt) antifcrro- 
magnetic correlation rather than the (tt/2, tt/2) antifcr- 
romagnetic correlation in bulk FeTe samples. Since the 
(tt, 0), (0, tt) antiferromagnetic scattering is natural from 
the band structure point of view, it is reasonable to start 
from the itinerant picture when studying the supercon- 
ductivity in FeTe films. 

We use the tight-binding model fitted to the band 
structure calculations for FeTe. The band structure has 
been partly confirmed by ARPES measurement 42 . Al- 
though it docs not capture the unusual magnetism in 
bulk FeTe, we hope it is an appropriate starting point to 
study the superconductivity in the FeTe thin films. Ex- 
cess iron exists in experimental samples of FeTe as in the 
case of FeSe, but will be ignored in our study. The Fermi 
surface, FRG flow and pairing form factor are presented 
in Fig. 5. As expected our results only show the (tt, 0) 
SDW correlation, not the (tt/2, tt/2) AFM in bulk FeTe. 

We note the following facts: (1) The Fermi surface 
shapes and their orbital contents(not shown) are very 
similar to that of LaFeAsO, except that the small hole 
pocket at T is much smaller, and has mainly xy orbital 
character. (2) The gap on electron pockets has substan- 



tial variations, but does not change sign. Like FeSe the 
large gap is concentrate in regions with dominant xy or- 
bital content. (3) The gaps on (tt, tt) hole pocket is large 
and has opposite sign from that on the electron pockets. 
Like FeSe, this suggests the root of pairing rests on the 
Cooper scattering between the (tt, tt) hole pocket and the 
electron pockets. (4) Interestingly, the gaps on the two 
central hole pockets are relatively small and have oppo- 
site sign. This suggests that these hole pockets play a 
secondary role in the superconducting pairing. The fact 
that the gap function reverse sign on the central hole 
pockets is a new feature that is absent in other iron based 
superconducting compounds. 

Real Space Representation of the Pairing Form 
Factors. Our FRG method can only calculate the pair- 
ing form factors around Fermi surfaces, and the pairing 
order parameter is obtained in band eigenbasis. To gain 
more intuitive picture, it is useful to have a real space 
picture of the pairing. 

Based on our experience and other theoretical works, 
we will focus on only the three t2 orbitals Xz,Yz,xy 
(here capital X, Y refer to the Fe-As directions, namely 
proper crystal axis) in the pairing order parameter. Off- 
site pairing between these three orbitals up to second 
neighbor on the Fe square lattice is included. According 
to the lattice symmetry the pairing order parameter has 
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FIG. 4: (Color online) Left panel: the Fermi surfaces of FeSe. Middle panel: The FRG flow of the interaction strength 
associated with the (n, 0) SDW and the superconducting pairing. Right panel: the pairing form factor. 
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FIG. 5: (Color online) Left panel: Fermi surfaces of FeTe. Middle panel: The FRG flow of the interaction strength associated 
with the (tt, 0) SDW and the superconducting pairing. Right panel: the pairing form factor determined by FRG. 



the following form, where *^ k)S = (icxz,-k,s icYz,-k, s Cxy,-k,a), super- 

script T means transpose, s =f, 4- labels the spin, i = 
I Ah (fez, k y ) A12 A 13 (k x ,k y ) \ \/^T, and the entries of the matrix are 

flkJ A12 A n (k y ,-k x ) A 13 {k y ,-k x ) * kjt 



-Ai 3 (k x ,k y ) -A ls (k y ,-k x ) A 



(3) 



An(k x , k y ) = A a cos(k x — k y ) + A v cos(k x + k y ) + A < - nn -'(cos k x + cos k y ), A12 = — Ai2(cos k x — cos k y ) 
A\z{k x , k y ) = A^' (sin k x — sinfcj,) + A° 3 sin^ — k y ), A33 = A < £ i n * > (cos k x + cosk y ) + A^ 1111 ' cosfe x cosfc, 



The pictorial representation of the eight fitting parame- 
ters A", A", A^, A 12 , A« Afc, A&\ and A 3 n 3 nn) 
arc illustrated in Fig. 6. 

Wc then project this order parameter onto the Fermi 
surfaces and compare with the FRG results. The fit- 
ting parameters are obtained by standard least square fit 
and listed in Tabic I. A similar real space representa- 
tion was recently obtained by Kariyado and Ogata from 
RPA results 43 . Note that although we have 24 points 
around one Fermi surface, due to symmetry of our 
model, there are actually only 24/8 = 3 independent data 
points to fit per Fermi surface. At this stage we have not 
been able to improve this due to limitation of computa- 



tional power. The fitting to LaFePO and FeSe are not 
very good, with large error estimates. This could come 
from our limited momentum space resolution, or that the 
pairing in these materials actually extend beyond second 
neighbor. So these results should be taken with caution. 

Nonetheless the differences between these materials 
are clearly visible, and are summarized below. (1) For 
LaFcAsO the intra-orbital pairing (A a and A< nn )), be- 
tween Xz-Xz or Yz-Yz, are strongest, therefore the 
gap is smallest in regions with dominant xy character 5 . 
(2) For LaFePO, the nearest- neighbor pairings, inter- 
orbital A12 between Xz-Yz and intra-orbital pairing 
A33 between xy-xy, are the strongest, which give a 
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FIG. 6: Pictorial representation of the fitting parameters and the three orbitals used. Empty and filled lobes in the orbital 
pictures indicate ± sign/symmetry of wave functions. Each of the fitting parameters, As, represents spin-singlet pairing between 
the two orbitals shown in the corresponding picture. 



nodal cos(k x ) + cos(k y ) gap. (3) In FeSe the gap is 
mainly from intra-orbital second neighbor xy-xy pair- 
ing Agg" 11 ^, which gives cos^) cos(k y ) type nodeless gap. 
The gap variation on electron pockets is mainly due to 
the variation of xy orbital content of the wave function. 
Thus it has the opposite anisotropy compared to our 
LaFeAsO result. (4) FeTe is very special, with dominant 

inter-orbital pairing (A12 and A^"-*) between Xz-Yz and 
X z/Yz-xy. The sign structure and variation of gap is a 
consequence of both this fact and the orbital content of 
the wave functions around Fermi surfaces. 

Conclusion We have so far studied the pairing of 
LaFeAsO, LaFePO, FeSe and FeTe (thin film) using the 
function renormalization group method. The most ro- 
bust picture that emerge from these studies is that pair- 
ing is due to the strong positive Cooper scattering be- 
tween the hole and electron Fermi surfaces. The positive 
value of such scattering is responsible for the tendency for 
the gap function to assume opposite sign on the electron 
and hole pockets. However, on each Fermi surface the 
average magnitude of the gap function, the degrees of its 
variation as a function of angle, and whether it has node 
or not depend on the details of the band wave function as 
well as the values of the local interaction in our model. 
The present result confirms our earlier picture 5 ' 20,31 of 
the antiferromagnctism-triggcred pair scattering being 
the mechanism of superconductivity in the iron-based 



compounds. We believe that such mechanism should 
also at work in systems at the brim of band-insulator to 
semi-metal phase transition (for example as a function 
of pressure). On the semi- metal side the charge count 
requires the electron and hole pockets to cover the same 
area (volume). If the centers of these pockets reside at 
time-reversal invariant fc-points, then intra-pocket pair- 
ing can take place. If the centers are not located at the 
time-reversal invariant points, then center-of-mass mo- 
mentum zero Cooper pair will requires inter-pocket pair- 
ing (time reversal ensure such k, —k pockets exists). It 
become increasingly clear with time that the highest T c 
compounds in the iron-based materials are not strongly 
correlating. This makes our method, the function renor- 
malization group an unbiassed, trustworthy method to 
study these systems. The neutron resonance and the 
quasiparticlc interference experiments raise our hope that 
perhaps after all these years we finally understand an un- 
conventional superconductor. 
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